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SUMMARY 


A  resonance  tube  method  has  been  used  to  measure  the  absorption  of  sound 
passed  through  a  column  of  propellant  gases  and  solid  particles  produced  by  the 
burning  of  small  samples  of  aluminized  propellant.  Aluminum  composition  was 
varied  to  observe  the  effect  of  the  increased  quantity  of  alumina  particles  on  sound 
absorption.  The  acoustical  absorption  coefficient  increased  with  increasing  aluminum 
concentration  in  the  propellant.  The  absorption  coefficient  increased  linearly  with 
the  frequency,  consistent  with  the  particulate  damping  theory  of  Epstein  and  Carhart.2 
This  theory  was  also  used  to  estimate  the  sizes  of  alumina  particles  produced  in  the 
experiment.  By  a  process  of  curve  fitting,  the  inferred  particle  radii  were  either 
0.80  or  3.  5  microns,  two  solutions  were  possible.  The  initial  average  particle 
radius  of  the  aluminum  was  2.  3  microns.  Despite  the  fact  that  a  statistical  distri¬ 
bution  of  particle  size  was  not  considered  in  cur ve -fitting  the  data,  the  reasonablenes 
of  the  calculated  particle  sizes  lends  further  support  to  the  Epstein  and  Carhart  theory. 
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NOMENCLATURE 
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velocity  of  sound 
frequency 

frequency  of  n  standing  waves 

resonant  frequency  range  at  half-power  of  peak  resonant  frequency 
tube  length 

order  of  resonance  or  integer  defining  harmonics,  cf.  Eq.  (4) 

number  density  of  particles 

effective  particle  radius 

absolute  temperature 

mole  fraction 
1/ 

(irf/v)  2  r,  dimensionless  particle  size 


Greek  Symbols 
6 
■n 


% 

G 

Subscript 

i 


ratio  of  gas  density  to  particle  density 
viscosity  of  gas 
kinematic  viscosity  of  gas 
average  molecular  weight 


acoustical  absorption  coefficient,  cm 

individual  components 
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INTRODUCTION 


Combustion  instability  in  solid  propellant  rockets  is  a  recurring  problem  which 
is  being  given  much  study.  As  has  become  clear,  it  is  a  complex  resonance  phenomenon 
intimately  connected  with  the  geometry  of  the  system,  the  nature  of  the  solid  pro¬ 
pellant,  and  a  variety  of  other  parameters.  Frequently,  it  has  been  found  that  the 
onset  of  this  condition  can  be  controlled  by  the  addition  of  certain  additives  to  the 
fuel.  The  mechanism  by  which  these  additives  work  is  by  no  means  certain:  they 
may  be  either  chemically  inert  or  may  actually  participate  in  the  combustion  process. 
For  instance,  one  common  additive --aluminum --forms  aluminum  oxide  which,  either 
as  solid  or  liquid  particles,  acts  to  suppress  the  instability,1  While  the  mechanism  by 
which  these  particles  act  is  not  clearly  established,  it  is  apparent  that  the  presence 
of  particles  in  the  acoustic  field  must  lead  to  the  absorption  of  field  energy  by  the 
particles  and  thereby  produce  a  damping  effect  on  the  undesirable  resonances.  If 
the  damping  effect  is  large  enough,  oscillations  will  not  occur. 

This  report  describes  a  series  of  experiments  undertaken  to  measure  the  con¬ 
tribution  of  the  alumina  particles  to  the  over -all  damping  in  a  solid  propellant  com¬ 
bustion  exhaust.  The  experiments  were  designed  to  separate  this  effect  from  the 
other  complicated  interactions  which  occur  in  a  real  motor. 


THEORETICAL  DISCUSSION 

In  a  rocket  motor  combustion  gas,  acoustical  energy  may  be  absorbed  by  at 
least  three  mechanisms:  absorption  by  walls,  gases,  and  particles  in  the  gas  phase. 

Acoustic  absorption  by  the  walls  of  a  tube  is  proportional  to  the  reciprocal  of  the 
tube  diameter  and  to  the  square  root  of  the  frequency.  It  is  also  relatively  small  and, 
as  the  experimental  values  show,  is  unable  to  explain  the  rather  large  damping 
obtained  here.  Absorption  by  the  gas  itself  is  usually  considered  to  be  due  to 
excitation  of  rotational  degrees  of  freedom  in  the  gas, therefore,  it  is  not  clear  how 
it  could  be  affected  by  addition  of  particles  such  as  alumina.  The  measured  values 

of  absorptivity  are  too  large  to  be  explained  by  this  phenomenon. 
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Epstein  and  Carhart2  have  developed  a  theory  for  the  attenuation  of  sound  by 
spherical  particles  suspended  in  a  gas.  In  this  theory,  both  the  viscous  and  thermal 
dissipative  losses  are  considered.  Since  the  thermal  losses  for  alumina  particles 
are  smaller  than  the  viscous  losses,  only  viscous  losses  will  be  considered  here. 

With  this  simplification  the  following  relationship  for  the  acoustic  absorption 
coefficient  is  obtained2: 

O’  =  (6irv/ c)  Nr  F(Z,  6)  (1) 

where: 

F(Z,6)  =  (1  +  Z)  16Z4  +  726  zi  +  8)g4  (1  +  2Z  +  2Z£) 

Z  E  ( trf  /  v)  ^ 1  r 

The  theoretical  expression  may  now  be  examined  to  understand  how  the  absorption 
coefficient  may  vary  with  particle  number  density,  N,  sound  frequency,  f,  and 
particle  radius,  r.  The  coefficient  O’  will  vary  linearly  with  the  particle  number 
density,  provided  (in  actual  combustion  systems)  the  particle  radii  do  not  differ 
appreciably  as  the  aluminum  concentration  is  varied.  In  order  for  N  to  vary 
linearly  with  the  aluminum  concentration,  it  is  assumed  that  the  aluminum  is 
burned  completely  to  its  oxide.  Sehgal3  has  reported  that  the  average  particle  sizes 
of  aluminum  particles  are  not  changed  appreciably  by  a  change  in  aluminum  content 
from  12  to  18  wt.%.  One  might  expect,  therefore,  that  O’  will  vary  linearly  with 
the  aluminum  concentration. 

By  assuming  a  constant  alumina  particle  radius  and  kinematic  viscosity  of  the 

propellant  gases,  it  may  be  shown  that  in  the  Z  range  of  interest  (Z  <  1)  O'*  f, 

where  the  exponent  n  varies  from  two  (Z<<  1)  to,  effectively,  zero  (Z  "large"). 

Therefore,  for  certain  fixed  particle  sizes,  it  may  be  possible  to  observe  an 

absorption  coefficient  which  varies  as  the  frequency  squared  to  one  which  is 

independent  of  sound  frequency.  In  particulars  (locally)  linear  variation  of  O’  with 

f  may  also  be  expected.  This  behavior  is  illustrated  in  Fig.  1  where  calculated 

values  of  F(Z)  for  the  particular  propellants  under  consideration  are  plotted  against 

Z  on  logarithmic  coordinates.  It  will  be  observed  that  the  lowest 

2 


(2) 

(3) 
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values  of  Z  give  a  line  with  a  slope  approaching  4  asymptotically,  while  for  higher 
values  of  Z,  F(Z)  approaches  zero  slope.*  Between  these  two  extremes,  the  curve 
shows  a  rather  sharp  break  near  Z  “  0.  01,  where  the  slope  changes  monotomically 
and  rapidly  from  almost  4  to  almost  0.  This  indicates  that,  for  fixed  values  of  r 
and  v,  O'  varies  as  f2  for  low  values  of  Z  but  becomes  essentially  independent  of 
frequency  at  somewhat  larger  values. 

The  Epstein -Car hart  formula  also  states  that  for  a  fixed  mass  fraction  of 
particles  the  absorption  coefficient  O'  varies  with  the  square  of  the  particle  size 
in  the  small  Z  extreme.  However,  for  larger  values  of  Z,  the  coefficient  O'  will 
vary  inversely  as  the  square  of  the  particle  size.  Hence  a  maximum  absorption 
coefficient  is  theoretically  predictable  for  an  optimum  particle  size.  For  alumina 
particles  suspended  in  the  propellant  product  gases  at  a  pressure  of  200  psi,  Horton 
and  McGie1  have  estimated  the  optimum  particle  radius  to  be  1.5  microns. 

The  parameters  N  and  6  are  sensitive  to  a  variation  in  temperature  because  ^ 

of  their  relation  to  gas  density.  When  all  variables  sensitive  to  temperature  change 
are  considered  (including  the  sound  speed  c)  the  absorption  coefficient  is  found  to 
vary  as  about  T_1.  + 

At  fixed  particle  mass  fraction  and  particle  size  pressure  does  not  affect  O' 
because  the  gas  particle  density  and  kinematic  viscosity  pres  sure -dependence 
cancel,  but  in  practice  pressure  may  nevertheless  influence  the  particulate  damping 
process.  Sehgal3  has  shown  that  the  pressure  level  in  an  experimental  strand  test 
motor  exerted  an  influence  on  the  average  particle  radii  of  alumina  produced  in 
these  burning  tests.  He  reported  that  the  average  alumina  particle  size  increased 
linearly  with  the  pressure  in  a  range  of  pressures  from  15  to  1000  psi.  Since 
sound  absorption  is  strongly  dependent  on  particle  size,  pressure  can  thus  exert 
a  strong  influence  on  the  particulate  damping  process  through  its  control  of  the  sizes 
of  particles  produed  in  the  burning  process.  One  is  led  to  conclude  that  the  mean 
pressure  in  the  rocket  motor  has  no  effect  on  the  gaseous  acoustical  damping  but 

* 

If  one  considered  still  larger  values  of  Z  then  F  would  ultimately  increase 
linearly  with  Z,  however,  these  values  of  Z  are  not  of  interest  here. 
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apparently  regulates  the  particulate  damping  process  via  the  above  mechanism.  Along 
with  this  interpretation  of  Sehgal's  work  goes  the  fact,  supported  by  the  Epstein- 
Carhart  theory,  that  there  is  an  optimum  particle  size  which  can  produce  a  maximum 
absorption  coefficient.  Consideration  of  this  would  mean  there  is  also  an  optimum 
pressure  which  is  associated  with  a  maximum  of  sound  absorption. 

In  this  discussion  of  particulate  damping,  it  has  been  shown  how  (T  might  be 
affected  by  aluminum  concentration,  frequency,  temperature,  pressure,  and 
particle  radius.  The  theory  predicts  a  linear  dependence  on  aluminum  concentration, 
a  varying  dependence  on  temperature  and  radius,  and  a  small  dependence  on  pressure 
when  pressure  does  not  influence  particle  size.  For  very  small  particle  sizes, 
acoustic  damping  can  become  small  compared  with  the  other  damping  processes. 
However,  within  the  range  of  particle  sizes  and  frequencies  experienced  in  this 
study  (0.  5  <  r  <  4  p;  0.  6  <  f  <  4  Kc/sec)  particulate  damping  can  be  the  dominating 
dissipative  mechanism. 

In  any  practical  case,  a  range  of  particle  sizes  will  exist,  so  that  this  theory-- 
which  is  for  single -sized  particles --will  require  modification. 


EXPERIMENTAL  TECHNIQUES 


It  is  usually  assumed  that  acoustic  instabilities  are  driven  by  complicated  non¬ 
linear  couplings  between  the  sound  field  and  the  flame  frcnt  in  the  rocket  motor.  As 
a  result,  it  is  most  difficult  to  separate  flame  front  effects  from  those  produced  by 
the  propellant  exhaust  gases  in  measurements  made  in  typical  motors.  In  order  to 
avoid  this  problem,  the  experiments  were  so  designed  that  the  flame  front  did  not 
contribute  to  the  acoustic  field.  Instead,  the  sound  waves  were  generated  by  an 
acoustic  driver,  while  the  solid  propellant  charge  served  only  to  supply  the  exhaust 
gases  whose  acoustic  absorption  it  was  desired  to  measure. 
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Open-Tube  Experiments 


A  simple  open-tube  technique  was  used  in  a  first  attempt  to  determine  the  relative 
sound  absorption  of  the  combustion  products  of  aluminized  and  nonaluminized  pro¬ 
pellants.  The  apparatus  consisted  of  a  brass  pipe  91.5  cm  long  with  a  probe  micro¬ 
phone  attached  midway  along  its  length.  A  60 -watt  acoustic  driver  was  coupled  to 
one  end  of  the  pipe,  while  a  replaceable  propellant  cartridge  was  located  near  the 
other  end  of  the  pipe  through  which  sound  waves  passed.  The  variation  in  sound 
intensity  with  time  as  measured  by  the  microphone  was  recorded  on  a  strip  chart 
as  combustion  gases  from  a  burning  solid  propellant  filled  the  tube.  These  tests 
were  run  at  a  number  of  frequencies,  not  necessarily  resonant  with  the  tube.  In  this 
way,  relative  sound  intensity  for  a  5  wt.  %  aluminum  concentration  were  determined 
to  be  15  db  below  those  of  a  nonaluminized  propellant.  With  higher  concentrations  of 
aluminum,  the  loss  rose  even  higher  but  appeared  to  level  off  at  10  wt.  %  aluminum, 
although  concentrations  as  high  as  20  wt.  %  were  tested. 

Accurate  interpretation  of  data  from  these  experiments  was  difficult,  however. 
Considerable  uncertainty  was  introduced  because  different  degrees  of  coupling 
between  driver  and  gas  and  between  gas  and  microphone  might  well  have  occurred  as 
the  gas  composition  changed.  The  physical  state  of  the  combustion  gases  was  not  well 
defined,  and  minor  resonances  and  background  noise  contributed  to  the  over -all 
uncertainty.  For  this  reason,  a  closed  resonance  tube  technique  was  used. 

Closed  Resonance  Tube  Experiments 

The  bulk  of  the  data  were  obtained  using  a  closed-tube  apparatus  adapted  from  an 
acoustic  experiment  of  Parker.4  This  technique  depends  upon  measuring  the  resonant 
frequency  for  standing  waves  in  a  tube  of  variable  length.  The  relative  width  of  a 
resonance  frequency  curve  is  measured,  which  is  equivalent  to  measuring  the  "Q" 
of  the  system.  Since  the  relative  width  of  the  frequency  curve  depends  primarily  on 
the  dissipation  of  energy  in  the  tube  and  only  weakly  on  other  parameters,  interpreta¬ 
tion  of  the  data  is  greatly  simplified. 
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A  rigid-wall  copper  tube  was  used  (Fig.  2)  90  cm  in  length  and  3.8  cm  in  diameter. 
A  1-gm  cube  of  solid  propellant  was  positioned  near  the  center  of  the  tube  and  ignited 
electrically.  At  one  end  of  the  tube,  a  short  section  of  narrower  tube  filled  with 
straight  steel  wires  served  as  a  coupling  tube  to  an  acoustic  driver.  This  coupling 
technique  afforded  a  uniform  distribution  of  sound  in  the  resonance  chamber  and  also 
provided  a  high  source  impedance.  *  The  length  of  the  tube  could  be  varied  to  obtain 
a  series  of  resonant  lengths  by  a  movable  piston  which  adjusted  the  resonance  fre¬ 
quency  of  the  assembly.  A  dynamic  crystal  microphone  was  mounted  on  this  piston 
and  mechanically  isolated  from  the  piston  by  a  cork  sleeve. 

Resonance  in  a  closed  tube  is  established  when  a  series  of  standing  waves  of 
frequency,  f  ,  are  contained  within  the  tube.  This  condition  occurs  when  the  fre¬ 
quency  is: 

f  =  nc/2L  (n  =  1,  2,  .  .  .)  (4) 

n 

When  the  propellant  starts  to  burn,  a  new  gaseous  medium  displaces  the  original  air 
and  a  new  series  of  standing  waves  appear  in  the  tube.  In  actual  experiments,  the 
resonant  frequency  shifts  to  somewhat  higher  values  because  of  the  lower  density  or 
the  higher  speed  of  sound  in  the  propellant  gases. 

Sound  absorption  was  measured  as  follows:  A  desirable  resonant  frequency  was 
first  established  in  air  and  the  driver  frequency  was  swept  automatically  at  a  rate 
of  16  cps  approximately  50  cps  on  either  side  of  the  resonant  frequency.  The  sweep 
was  from  low  to  high  frequency  and  then  reversed  to  remove  any  hysteresis  problems. 
The  microphone  response  was  recorded  with  time  and  displayed  a  peak  output  each 
time  the  frequency  passed  through  resonance.  Determination  of  the  acoustical 
absorption  coefficient  involved  measurement  of  peak  width  at  the  "half -power"  points 
of  the  resonant  curve  together  with  the  determination  of  the  resonant  frequency. 

The  acoustical  absorption  coefficient  is  then  given  by  the  expression: 

* 

Other  forms  of  coupling,  including  direct  coupling,  were  also  tried  with  little  real 
difference  in  results. 
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(T  =  ttA  f/ c 

or  since  c  =  2f  L/n 
n 

O'  =  ( irn/ 2 L)  (Af/f^)  (5) 

where  Af  is  the  resonant  frequency  range  at  the  half -power  points.  The  order  of 
resonance,  n,  was  determined  from  the  resonance  peak  measurements  in  air,  since 
f^,  c,  and  L  are  then  known  in  Eq.  (4).  Knowing  n,  the  value  of  c  for  any  particular 
gas  or  gas  mixtures  was  determined. 

A  block  diagram  of  the  instrumentation  is  given  in  Fig.  3.  Driver  input  was 
obtained  through  a  power  amplifier  from  the  signal  generator  of  a  Hewlett-Packard 
Model  500  (1%  of  full  scale,  lOOKc)  wave  analyzer  whose  frequency  was  recorded 
on  one  channel  of  a  dual -channel  Sanborn  Model  322  recorder.  Output  of  the  micro¬ 
phone  passed  through  a  wave  analyzer  (acting  as  a  filter)  and  was  recorded  on  the 
other  channel  of  the  recorder.  The  recordings  therefore  gave  simultaneous  measure¬ 
ments  of  frequency  and  microphone  response  at  all  times. 

Upon  ignition  of  the  solid  propellant,  the  resonant  frequency  of  the  tube  could  be 
observed  to  change,  as  well  as  the  relative  peak  shape  and  peak  heights.  A  series 
of  such  peaks  were  recorded  for  each  run  to  determine  the  time  response  of  the  system. 
Propellant  samples  tested  were  composite  mixtures  of  ammonium  perchlorate 
oxidizer  and  polystyrene  fuel.  In  order  to  incorporate  varying  amounts  of  aluminum, 
the  basic  composition  (78  wt.  %  NH4C104-Z2  wt.  %  polystyrene)  was  changed  in  the 
manner  shown  in  Table  I.  This  composition  change  produced  a  310°K  change  in  the 
calculated  adiabatic  flame  temperatures  at  the  two  extremes  (0  and  10%  aluminum). 


EXPERIMENTAL  RESULTS 

Acoustical  absorption  coefficients  were  evaluated  by  a  series  of  runs  in  the 

closed-tube  experiments.  The  first  experiments  were  done  with  the  nonaluminized 

propellant  and  are  presented  in  Fig.  4  with  results  from  Horton  and  McGie1  and 

theoretical  curves.  The  Horton  and  McGie  data  were  obtained  at  14  atm  total 

pressure,  and  our  data  were  obtained  at  1  atm.  The  propellant  employed  by  Horton 

and  McGie  was:  80  wt.  %  ammonium  perchlorate,  18  wt.  %  polybutadiene  polyacrylic 
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acid,  and  2  wt.  %  copper  chromite.  The  theoretical  curves  were  computed  with  the 
Helmholtz -Kirchoff  equation  for  a  plane  wave  in  a  tube  with  damping  using  the 
calculated  adiabatic  flame  temperatures  and  corresponding  gas  compositions. 

The  two  experimental  results  are  rather  similar.  The  deviation  from  theory 
is  not  understood  but  is  probably  due  in  part  to  a  lower -than-calculated  temperature. 

In  order  to  further  understand  the  interrelationship  of  the  various  parameters 
involved  in  these  measurements,  a  group  of  studies  was  conducted  to  determine 
the  change  in  absorption  with  time  in  the  resonance  tube.  Studies  for  5  wt.  % 
aluminum  samples  in  both  an  open  and  a  closed  tube  and  0  wt.  %  aluminum  in  an 
open  tube,  are  presented  in  Fig.  5;  in  the  closed  tube,  pressures  rose  to  2.  04  atm 
following  propellant  ignition,  whereas  in  the  open  tube  the  pressure  remained  at 
atmospheric.  It  was  observed  that  an  initial  sharp  drop  in  absorption  occurred  in 
the  open  tube  well  within  the  first  100  sec  followed  by  a  much  slower  decrease  over 
an  extended  period.  The  sharp  drop  was  not  nearly  so  evident  in  the  closed  tube, 
with  the  final  portions  of  the  two  curves  having  essentially  identical  slopes.  No 
explanation  for  the  difference  in  the  initial  portion  of  the  curves  is  proposed. 
Experimental  details  make  measurement  of  Af  difficult  at  any  time  earlier  than 
10  sec  after  burning  ceases,  so  that  the  initial  portion  of  the  curve  is  somewhat 
dubious.  The  central  resonant  frequency,  3100  cps,  did  not  change  appreciably 
over  the  indicated  time  period  so  that  the  time  variation  of  Af  corresponds  to  a 
time  variation  of  the  absorption  coefficient  d  .  In  this  time,  the  Af  value  would 
have  shifted  about  7  cps.  This  could  represent  an  error  of  about  10%  in  the  measure¬ 
ments  reported  here  and  would  always  represent  a  lower  absorption  coefficient  than 
the  true  value.  It  does  appear  that  in  the  open-tube  experiments  the  absorption 
coefficient  decreased  in  an  exponential  manner  for  the  first  60  sec  with  a  half-life 
of  about  45  sec.  For  the  next  15  min  the  absorption  decreased  exponentially  with 
a  half-life  of  about  700  sec.  The  only  effect  of  burning  at  a  slightly  higher  pressure 
was  to  make  the  initial  period  of  drop  of  doubtful  significance,  without  much  effect 
on  the  long-time  rate  of  change.  It  is  not  clear  whether  the  initial  drop  occurred 
more  rapidly  or  not  at  all. 
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The  change  with  time  may  be  due  to  several  possible  effects:  (1)  a  decrease  in 
the  total  number  of  particles  interacting  due  to  settling  or  agglomeration,  (2)  a 
change  in  average  particle  radius  due  to  a  differential  settling  rate  among  the 
various  sizes  of  particles  present,  (3)  a  change  in  average  radius  of  particles 
due  to  agglomeration  of  particles,  and  (4)  a  change  in  absorption  due  to  changes 
in  ambient  conditions  of  the  exhaust  gas,  i.  e.  ,  changes  in  temperature,  pressure, 
composition,  etc.  (Note  in  Fig.  5  that  the  rate  of  decay  is  the  same  with  and  without 
aluminum. ) 

While  it  will  be  shown  that  CT  will  be  most  sensitive  to  changes  in  particle  radius, 
very  large  changes  in  temperature  will  still  occur  as  the  particles  cool  off,  so  that 
temperature  could  also  be  significant.  It  is  evident  that  the  change  in  time  of  the 
absorption  coefficient  will  be  a  major  complication  in  this  experiment.  Should 
future  work  be  done  with  this  technique,  this  complication  will  have  to  be  given 
more  attention  than  was  possible  here. 

4 

The  results  with  various  concentrations  of  aluminum  are  given  in  Table  II. 

The  absorption  coefficients  are  plotted  against  aluminum  concentration  in  Fig.  6. 

The  values  of  Af  and  f  represent  an  average  of  several  experimental  runs.  The  t 

mean  deviations  given  with  Af,  f  and  O’  increases  significantly  with  increasing 
aluminum  content,  particularly  at  higher  resonant  frequencies. 

The  absorption  of  sound  is  expected  to  increase  linearly  with  aluminum  con¬ 
centration  if  the  particle  size  remains  constant,  i.  e.  ,  the  number  density  of 
particles  varies  linearly  with  aluminum  concentration.  It  has  been  observed  by 
H.  Silla  of  this  Laboratory  that  aluminized  propellants  burn  inefficiently  in  1  atm 
of  nitrogen.  Rough  estimates  indicate  that  only  60  wt.  %  of  the  aluminum  is  burned. 

At  high  (about  10%)  concentrations  of  aluminum,  the  efficiency  is  less  than  at 
lower  concentrations.  Thus,  the  lack  of  linearity  in  Fig.  5  may  be  due  to  inefficient 
combustion  of  aluminum  which  would  be  expected  to  produce  a  different  particle  size. 
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COMPARISON  OF  EXPERIMENT  AND  THEORY 

Equation  (1)  was  used  to  calculate  a  theoretical  O'  in  order  to  compare  with 
measured  values  from  the  resonance  tube  technique. 

The  various  properties  of  the  propellant  gases  (density,  viscosity,  particle 
number  density,  etc.)  were  derived  from  thermodynamic  calculations  of  the 
adiabatic  flame  temperature.  The  viscosity  of  the  gaseous  mixtures  was  calculated 
from  the  approximate  formula 


r]  =  (lA  ) 


X.  %.  Tl. 

I  l  1 


(6) 


The  individual  gas  viscosities,  rj^,  were  taken  from  the  viscosity -temperature  tables 
assembled  by  Svehla.  5  The  alumina  particle  density  was  assumed  to  be  4.  0  gm/ cm3, 
the  same  as  the  bulk  density  of  A1203. 


Table  III  lists  the  calculated  gas  properties.  Expressions  for  N  and  O'  are  also 
given.  The  expression  for  O'  is  simplified  in  terms  of  the  function  F(Z)  and  r,  both 
unknown  quantities.  The  expression  for  O'  is: 

CT  =  constant  •  F(Z)/r2  (7) 

The  constant  in  Eq.  (7)  is  the  term  (6ttv/c)N  appearing  in  Eq.  (1),  which  has  now 
been  calculated. 

The  function  of  F(Z)  may  be  estimated  by  assigning  arbitrarily  chosen  values 
of  Z  and  the  calculated  ratio  of  gas  to  particle  density,  6  .  Figure  1  describes  the 
variation  of  the  function  F(Z)  in  terms  of  Z.  Note  the  curve  has  a  steep  slope  where 
Z  has  values  of  5  x  10-3  and  lower.  As  Z  increases,  the  slope  of  the  curve  pro- 

V 

gressively  decreases  to  values  near  zero.  Since  Z  depends  linearly  on  r  and  f  , 

O'  may  be  expected  to  vary  with  frequency  when  r  is  constant,  according  to  the  value 
of  the  slope  of  the  F(Z)  vs.  Z  curve.  If  the  slope  is  4,  C  will  vary  as  the  frequency 
squared.  A  slope  of  2  will  produce  a  linear  variation  with  frequency.  A  slope  of 
zero  means  that  O'  is  independent  of  the  frequency. 


A 
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The  observed  absorption  coefficients  are  plotted  as  O’  vs.  f  in  Figs.  7  to  10. 

The  observed  slopes  of  unity  displayed  by  the  experimental  curves  show  that  the 
composition  of  the  gases  and  particles  is  such  as  to  produce  a  2  with  a  value 
nearly  10-2,  where  0  varies  as  Z2  or  f  (Fig.  1). 

It  was  possible  to  calculate  O’  using  Eq.  (l)  and  Fig.  1.  In  this  calculation, 
various  values  of  particle  radii  were  assumed  and  the  experimental  frequencies 
were  used.  One  experimental  O'  was  then  fitted  (at  a  frequency  of  3070  cps)  with 
a  calculated  O' .  Once  this  experimental  point  was  fitted,  r  was  fixed  and  other  O’ 's 
were  calculated  corresponding  to  other  frequencies,  and  these  absorption  coefficients 
were  compared  with  those  actually  measured  by  the  resonance  tube  technique. 

Figures  7  and  8  show  this  comparison  of  observed  and  calculated  0  as  a  function 
of  frequency;  both  variables  are  plotted  on  logarithmic  coordinates.  It  is  apparent 
that  the  experimental  points  form  a  straight  line  of  slope  unity,  which  passes 
generally  between  the  two  curves  calculated  from  the  theoretical  expression  for  d  . 
Since  the  two  calculated  curves  represent  assumed  small  and  large  particle  sizes 
(both  of  which  fit  the  experimental  CF  at  3070  cycles/sec),  it  is  impossible  to  decide 
on  which  particle  size  corresponds  to  reality. 

Figures  7  and  8  display  the  calculated  values  of  CT  for  propellant  gases  assumed 
to  be  near  the  adiabatic  flame  temperatures,  2730°K  and  2870°K,  for  the  5  and  10 
wt.  %  aluminum  concentrations,  respectively.  This  temperature  may  be  too  high, 
since  appreciable  cooling  of  the  propellant  gases  must  have  taken  place  in  the 
resonance  tube. 

To  observe  the  effect  of  a  lower  temperature,  the  above  calculation  was  repeated 
for  a  temperature  of  1000°K,  where  corrections  were  applied  to  viscosity,  gas 
density,  velocity  of  sound,  and  all  other  parameters  which  might  be  affected  by  a 
decrease  in  temperature.  Figures  9  and  10  show  the  results  of  this  calculation. 

There  is  somewhat  closer  agreement  of  the  calculated  and  experimental  curves 
for  the  two  propellants,  but,  again,  the  experimental  curve  is  approximately  an 
average  of  the  curves  calculated  from  theory. 
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It  is  again  impossible  to  determine  which  particle  size  represents  reality.  A 
calculation  was  also  carried  out  at  500°K  with  no  improvement  in  the  agreement  of 
O’  (observed)  and  CT  (calculated). 

For  the  5  wt.  %  aluminum  concentration,  the  small  particle  radius  was  calculated 
to  be  0.  95  micron  and  the  large  particle  radius  was  calculated  to  be  3.  4  microns. 

The  10  wt.  %  aluminum  concentration  produces  a  greater  difference  in  particle  radii, 

0.  55  and  3.  8  microns. 

Particle  number  densities  vary  from  104  to  107  particles/cm3  for  the  large  and 
small  particles,  respectively.  Horton  and  McGie3  have  calculated  a  particle  number 
density  of  107  cm-3  for  a  propellant  containing  1%  Al,  where  r  was  0.  6  microns.  A 
comparison  with  their  calculation  of  particle  number  density  is  offered  only  to  show 
that  the  calculations  given  here  are  reasonable.  As  noted  earlier  larger  particle 
number  densities  obtained  for  smaller  particle  radii  cannot  be  interpreted  to  mean 
that  a  larger  absorption  of  sound  is  theoretically  possible  for  smaller  particles. 

Size  must  also  be  considered  in  the  light  of  the  previously  discussed  existence  of  an 
optimum  particle  size  producing  a  maximum  of  sound  absorption.  This  optimum 
radius,  as  calculated  by  Horton  and  McGie, 1  is  1.5  microns  for  alumina  particles, 
for  propellant  compositions  of  unusually  low  concentrations  of  aluminum  (0  to  1.  5  wt.%) 
and  for  combustion  carried  out  at  a  pressure  of  200  psi.  In  our  investigation,  in 
which  propellants  contained  as  much  as  10  wt.%  aluminum  and  were  burned  at  1  atm 
pressure,  one  might  reasonably  expect  a  larger  particle  size  to  be  produced  and 
also  that  the  optimum  particle  radius  would  be  somewhat  larger  than  1.  5  microns. 

The  two  possible  average  particle  radii  inferred  from  these  experiments  are  0.80 

and  3.  5  microns  averaged  from  the  results  reported  in  Figs.  7  through  10: 

Small  radii  Large  radii 

1.  7  microns  3.  0  microns 

1.0  4.0 

0.95  3.4 

0.55  3.8 

The  initial  average  particle  radius  of  the  aluminum  used  was  2.  3  microns.  These  fall 
well  within  the  domain  of  particle  sizes  which  have  been  observed  directly3  and  this 
in  itself  gives  confidence  to  the  experiment  and  procedure. 

The  experimental  work  should  be  repeated  more  carefully  and  analyzed  using  a 
theory  valid  for  a  statistical  size  distribution  of  particles.  These  were  only  pre¬ 
liminary  experiments  which  demonstrate  the  possible  utility  of  the  method. 
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TABLE  I 


PROPELLANT  COMPOSITIONS  AND  CALCULATED  FLAME  TEMPERATURE 


Al,  wt.% 

NH4C104,  wt.% 

P-13,  *wt.% 

T,  °K 

Gas  visco 
poise 

sity  Gas  density, 

1 0  x  gm/  cm 

0.  0 

78.  0 

22.  0 

2560 

8.  09 

1.  12 

2.  0 

76.  4 

21.6 

— 

— 

— 

4.  0 

74.  9 

21.  1 

____ 

— 

— 

5.  0 

74.  1 

20.9 

2730 

8.63 

1.  09 

8.  0 

71. 7 

20.  3 

— 

— 

— 

10.  0 

70.  2 

19.8 

2870 

6.  36 

1.  07 

* 

Polystyrene 
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TABLE  II 


EFFECT  OF  ALUMINUM  CONCENTRATION  ON  ABSORPTION  COEFFICIENT 


Aluminum, 

Frequency, 

Frequency  width, 

Velocity  of  sound. 

CT  ,  abs  coeff, 

wt.  % 

cps* 

cps* 

104  cm/ sec* 

1C3  cm-1* 

0.  0 

645  ±  7 

26. 0  ±  3.  0 

3.  59  ±  0.  04 

2.  26  ±  0.  32 

2.  0 

643  ±  10 

26. 0  ±  3.  0 

3.  60  ±  0.  05 

2.  24  ±  0.  32 

4.  0 

650  ±  4 

29.  0  ±  7.  0 

3.  64  ±  0.  02 

2.48  ±  0.61 

5.  0 

649  ±  2 

30.  0  ±  8.  0 

3.63  ±  0.  01 

2.  58  ±  0.  65 

8.  0 

644  ±  0 

29. 0  ±  3.  0 

3.  60  ±  0.  00 

2.  56  ±  0.  25 

10.  0 

658  ±  4 

33.  0  ±  5.  0 

3.68  ±  0.  02 

2.85  ±  0.  40 

0.  0 

1570  ± 10 

23.  0  ±  2.  0 

3.  53  ±  0.  01 

2.  00  ±  0.  74 

5.  0 

1560  ±  20 

63.  0  ±  0.  5 

3.  48  ±  0.  02 

5.  53  ±  0.82 

10.0 

1560  ±  0 

65.  1  ±  0.  0 

3.  50  ±  0.  0 

5.  70  ±  0.  29 

0.  0 

3077  ±  21 

44.  0  ±  5.  0 

3.  44  ±  0.  02 

4.  17  ±  0.  41 

2.  0 

3090  ±  7 

50.  0  ±  3.  0 

3.  45  ±  0.  01 

4.  51  ±  0.  28 

5.0 

3050  ±  0 

83.  0  ±  2.  0 

3.  42  ±  0.  00 

7.  61  ±  0.  16 

o 

• 

CO 

3070  ±  0 

105  ±7.0 

3.  44  ±  0.  00 

9.  64  ±  0.  60 

10.  0 

3060  ±  10 

107  ±7.0 

3.  43  ±  0.  01 

9.  70  ±  0.  68 

0.  0 

4020  ± 10 

79.  3  ±  1.  1 

3.68  ±  0.  02 

6.68  ±  0.  20 

5.  0 

4035  ±  35 

15.  5  ±  2.  0 

3.  72  ±  0.  02 

12.9  ±0.30 

(  ±  mean  deviation) 
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TABLE  III 

PROPERTIES  OF  PROPELLANT  GASES  USED  IN  THE 


CALCULATION  OF  C 

Composition  of  Propellant 

74.  1  (AP)  -20.  9 

70.  2  (AP)  -9.  8 

sytrene  fuel 

styrene  fuel 

5  A1 

10  A1 

Flame  Temperature,  T 

2730°K 

2870°K 

Gas  Density,  gm/ cm3,  p 

1.  09 

1.  07 

Viscosity,  (poise)  r| 

8.63 

6.  36 

Kinematic  Viscosity,  cm2/sec,  v 

7.9 

5.94 

Speed  of  Sound,  cm/sec,  c 

1.  07  x  10s 

1.  08  x  105 

Ratio  of  gas  to  Particle  Density, 

2.  73  x  lO"6 

2.68  x  lO'5 

Mass  of  Aluminum  Oxide,  gm,  ms 

0.  095 

0.  190 

Mass  of  Gases,  gm,  mg 

0.905 

0.810 

Particle  Number  Density  cm'3,  N 

6.  83  x  lO'7 

1.49  x  lO"6 

r3 

r3 

Acoustic  Absorption  Coefficient,  calc. 

9.45  x  lO"10  F(Z) 

1.  58  x  lO"9  F(Z) 

rz 

z 

r 

* 
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2  4  6  8  |0-2  2  4  6  8  |QH 

Z 

FIG.  1  VARIATION  OF  f(Z)  WITH  THE  DIMENSIONLESS  DISTANCE 
IN  EPSTEIN-CARHART  FORMULA 

Assumed  t  wt.  %  Al,  T  =  2700°K,  P  =  1  atm, 
corresponding  to  6  =  2.  73  x  10"5 
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FIG.  ACOFSiK  RFSOXANCI 


FIG.  3  FLOW  CHAR  T  OF  ACOUSTICAL  AND  ELECTRONIC  EQUIPMENT 


FREQUENCY,  kilocycles/sec 
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FIG.  5  VARIA  TION  IN  ABSORPTION  COEFFICIENT  WITH  TIME 

AFTER  IGNITION 
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FREQUENCY  ,  cycles/sec. 

FIG.  7  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  ACOUSTIC 

ABSORPTION  COEFFICIENTS 

Theoretical  curves  based  on  particle  radii  to  fit  experiment  at  3070  cps 

5  wt.  %  Al,  T  =  2730°K 
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ACOUSTIC  ABSORPTION  COEFFICIENT,  CT,  cm” 


6S-60 


4  6  8  |08  2  4  6  8  |04 

FREQUENCY  ,  cycles/sec. 


FIG.  8  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  ACOUSTIC 

ABSORPTION  COEFFICIENTS 

Theoretical  curves  based  on  particle  radii  to  fit  experiment  at  3070  cps 

5  wt.  %  Al,  T  =  2730°K 
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4  6  8  j03  2  4  6  8  |q4 

FREQUENCY,  cycles/sec. 

FIG.  9  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  ACOUSTIC 

ABSORPTION  COEFFICIENTS 

Theoretical  curves  based  on  particle  radii  to  fit  experiment  at  3070  cps 

5  wt.  %  Al,  T  =  1 00U°K 

2  5 


Theoretical  curves  based  on  particle  radii  to  fit  experiment  at  3070  cps 

10  ,wt.  %  Al,  T  =  1  000°K 
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1  ABSTNAC1' 

A  resonance  tube  method  has  been  used  to  measure  the  absorption  of  sound 
(passed  through  a  column  of  propellant  gases  and  solid  particles  produced  by  the 
Jburning  of  small  samples  of  aluminized  propellant.  Aluminum 

q'.aiiti',  of  alumina,  particles  on  sound 
bsorption.  The  acoustical  absorption  coefficient  increased  with  increasing 
luminum  concentration  in  the  propellant.  The  absorption  coefficient  increased 
inearly  with  the  frequency,  consistent  with  the  particulate  damping  theory  of 
pstein  and  Carhart.  This  theory  was  also  used  to  estimate  the  sizes  of  alumina 
barticles  produced  in  the  experiment.  By  a  process  of  curve  fitting,  the  inferred 
particle  radii  were  either  0.80  or  3.  5  microns,  two  solutions  were  possible.  The 
nnitial  average  particle  radius  of  the  aluminum  was  2.  3.  Despite  the  fact  that  a 
statistical  distribution  of  particle  size  was  not  considered  in  curve -fitting  the  data, 
(the  reasonableness  of  the  calculated  particle  sizes  lend  further  support  to  the 
(Epstein  and  Carhart  theory. 
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